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^ I Abstract: 

Applying the Hamilton-Jacobi method beyond the semiclassical approximation prescribed in 
. [12j for the scalar particle. Hawking radiation as tunneling of Dirac particle through an event 

I horizon is analysed. We show that, as before, all quantum corrections in the single particle ac- 
tion are proportional to the usual semiclassical contribution. We also compute the modifications 
to the Hawking temperature and Bekenstein-Hawking entropy for the Schwarzschild black hole. 
Finally, the coefficient of the logarithmic correction to entropy is shown to be related with the 
QJ^' trace anomaly. 



Oh: 



I Semiclassical methods of modeling Hawking radiation as a tunneling effect were developed 

■ over the past decade and have generated a lot of interest [Il[2l[3llll[5l[6l[7l[8l[9l [lO] . Prom 
, this approach an alternative (intuitive) way of understanding black hole radiation emerged. 
However, most of the calculations in the literature pi El HJ [5l [U [71 [U EJ [TO] have been performed 
just for scalar particles. Since a black hole can radiate all types of particles like a black body, 
. the emission spectrum should contain particles of all spins. Therefore a detailed study of spin 
00 ! one-half particle emission is necessary. Although there exist some computations [11] in this 
' context, these are confined to the semiclassical approximation and do not consider quantum 
corrections. 

^ , In our previous work [12], we formulated the Hamilton-Jacobi method of tunneling beyond 

?H I semiclassical approximation by considering all the terms in the expansion of the one particle 
action for a scalar particle. We showed that the higher order terms are proportional to the 
semiclassical contribution. This result, together with properties of conformal transformations, 
eventually led to corrected expressions for thermodynamic variables of a black hole. It is not 
obvious whether a similar analysis is valid for the case of spin-half fermion tunneling. This issue 
is addressed here. 

In this paper we will discuss the Dirac particle tunneling beyond semiclassical approxima- 
tion employing the Hamilton-Jacobi method suggested in [12]. We will explicitly show that 
the higher order terms in the single particle action are again proportional to the semiclassical 
contribution. By dimensional argument the form of these proportionality constants, upto some 
dimensionless parameters, are determined. In particular for Scwarzschild spacetime, these are 
given by the inverse powers of the square of the mass of the black hole, because in this case, the 
only macroscopic parameter is mass. Using the principle of "detailed balance" [21 [9] the modified 
Hawking temperature is identified. Then the corrections to the Bekenstein-Hawking area law 
are derived by using the Gibbs form of first law of thermodynamics. Interestingly, the leading 
order correction to the entropy is the logarithmic of the semiclassical entropy which was found 
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earlier in[I31[Tl[I51[ISl[I71[ISl[ia[2ni[2IlE2]. Finally, using a constant scale transformation to 
the metric, we show that the coefficient of the logarithmic correction is related to trace anomaly. 



Our method involves calculating the imaginary part of the action for the (classically forbid- 
den) process of s-wave emission across the horizon which in turn is related to the Boltzmann 
factor for emission at the Hawking temperature. We consider a massless Dirac particle in a 
general class of static, spherically symmetric spacetime of the form 



9{r) 



(1) 



where the horizon r = rj^ is given by /{rn) = di^n) = 0. The massless Dirac equation is given 
by 

n^V.V = (2) 



where for this case the 7 matrices are defined as. 
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The covariant derivative is given by. 
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rsinO \ cr 



(3) 



7a, 7/3 



(4) 



Since for radial trajectories only the (r — t) sector of the metric ([T|) is important, ([2]) can be 
expressed as 



Here the required nonvanishing connections are 



^ tt 



111 rt =IL 
2 ' 2f 



Therefore, under the metric ([T]), the Dirac equation ([2]) reduces to 
and the matrix form of E^t from ([3D and (jH) is given by 



-Tt 
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To solve d?]) we employ the following ansatz for the spin up (i.e. +ve r-direction) and spin down 
(i.e. -ve r-direction) ijj as 



V't(t,r) 



/ A(t,r) \ 


B{t,r) 

V ; 



exp 



(9) 



V'i(t,r) 



/ \ 

C{t,r) 


V D{t,r) J 



exp 



(10) 



where I{r, t) is the one particle action which will be expanded in powers of h. Here we will only 
solve the spin up case explicitly since the spin down case is fully analogous. Substituting the 
ansatz ([9]) in ([7]), we obtain the following two equations: 



(11) 



Vf 



0. 



(12) 



Since the last terms within the first bracket of the above equations do not involve the single par- 
ticle action, they will not contribute to the thermodynamic entities of the black hole. Therefore 
we will drop these two terms. Now taking 1-^ = 1 and expanding /, A and B in powers of h, we 
find, 



Iir,t)=Io{r,t)+J2fi'liir,t) 

i 

A = AQ + Y,^"Ai- B = Bo + Y.^'^^- 



(13) 



where i = 1,2, 3, In these expansions the terms from 0{h) onwards are treated as quantum 

corrections over the semiclassical value /q, Aq and Bq respectively. Substituting (fTB|) in (fTT|) and 
(jl2p and then equating the different powers of h on both sides, we obtain the following two sets 
of equations: 



Set I : 



Vf 



AodJo + ^Bodrlo = 



-jjAodth + -jjAidth + VgBodrh + ^Bidrh = 

-^Aodth + -^Aidth + -^A2dth 
y/gBodrh + y/gBidrh + ^/gB2drh = 



(14) 
(15) 

(16) 



and so on. 
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Set II : ^1° : 



-^Bodth - -^BidJo + VgAodrh + ^Aidrh = 

-^Bodth - -^B^^th - -^B2dth 
+ y/gAodrh + ^/gAidrh + ^Aidrh = 



(17) 
(18) 

(19) 



and so on. 

Equations (jl4p and (jl7p are collectively known as the semiclassical Hamilton-Jacobi equations 
for a Dirac particle. Since the metric ([1]) is stationary it has timelike Killing vectors. Thus we 
will look for solutions of ()14p and (|17p which behave as 

Io = wt + T^(r), (20) 

where uj is the energy of the particle. Substituting this in and pT|) we obtain, 

'-^u, + B^^W'{t) = Q 

-'^u: + Ao^W'{r) = ^. (21) 

These two equations have two possible solutions: 

dr 



Ao = iBo; W+{r)=uj f 

Jo 



\/f{r)g{r) 



Ao = -iBo; W.{r) = -u; (22) 

where W+(VF_) corresponds to ingoing (outgoing) solutions. The limits of the integration are 
chosen such that the particle goes through the horizon r = vh- Therefore the solution for Io{r, t) 
is 

Io{r,t)=ut±u / (23) 
-'o \/f{r)g{r) 



Now, it is interesting to note that using (p2|) and (p3|) in the equations of Set I and Set II 
simultaneously and then solving we get relations connecting different orders in the expansion of 
A with those of B: 

Aa = ±iBa (24) 

where a = 0, 1, 2, 3, These lead to a simplified form of all the equations in Set I and Set II 

as. 



dtia = ±VfgdrIa (25) 

i.e. the functional form of the above individual linear differential equations is same and is iden- 
tical to the usual semiclassical Hamilton-Jacobi equations (|14p and (jl7p . Therefore the solutions 
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of these equations are not independent and /j's are proportional to /q- A similar situation hap- 
pened for scalar particle tunneling [12]. Since Iq has the dimension of H the proportionality 
constants should have the dimension of inverse of Again in the units G = c = = 1 the 
Planck constant h is of the order of square of the Planck Mass Mp and so from dimensional 
analysis the proportionality constants have the dimension of M~^* where M is the mass of black 
hole. Specifically, for Schwarzschild type black holes having mass as the only macroscopic pa- 
rameter, these considerations show that the most general expression for I, following from (113p . 
valid for (p5]) . is given by, 



/(r,t) = (l+^A-^)/o(r,t). 



(26) 



where /3j's are dimensionless constant parameters. 

The above analysis shows that to obtain a solution for /(r, t) it is therefore enough to solve for 
Iq (r, t) which has the solution of the form (|23p . In fact the standard Hamilton- Jacobi solution 
determined by this Ioir,t) is just modified by a prefactor to yield the complete solution for 
I{r,t). Substituting in we obtain 



KrA) 



1 + 



M2« 



dr 



\/f(r)g(r) 



Therefore the ingoing and outgoing solutions of the Dirac equation 
metric ([T|) is given by exploiting ([9]) and (l2j 



m 



exp 



dr 



and 



^out ~ exp 



■)( 



uot — UJ 



\/ f{r)g{r) 



dr 



^(1 + EA^2.^V"^ ~7o yWW) 



(27) 

under the background 
(28) 

(29) 



Now for the tunneling of a particle across the horizon the nature of the coordinates change. The 
sign of the metric coefficients in the (r — t) sector is altered. This indicates that 't' coordinate 
has an imaginary part for the crossing of the horizon of the black hole and correspondingly there 
will be a temporal contribution to the probabilities for the ingoing and outgoing particles. This 
has similarity with [23] where they show for the Schwarzschild metric that two patches across 
the horizon are connected by a discrete imaginary amount of time. 

The ingoing and outgoing probabilities of the particle are, therefore, given by. 



m 



ml 



exp 



dr 



\/f{r)g{r] 



and 



Pr 



out 



exp 



2 K 

-(1 + ^ A^)(u;Im t - wim 



dr 



Vf(r}gir] 



(30) 



(31) 



Now the ingoing probability Pjj^ has to be unity in the classical limit (i.e. h ^ 0) - when there 
is no reflection and everything is absorbed - instead of zero or infinity jl2j.Thus, in the classical 
limit, (1301) leads to. 



Im t 



-Im 



dr 



\/f{r)g{r) 



(32) 
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From the above one can easily show that Im t = —2ttM for the Schwarzschild spacetime which 
is precisely the imaginary part of the transformation t ^ t — 2i'KM when one connects the two 
regions across the horizon as shown in [23J . Therefore the probability of the outgoing particle is 

^out~exp[--c.(l + EA^jlmy^ 77(^1^ ^ ' ^''^ 
Now using the principle of "detailed balance" O [9] 

Pout = exp( - = exp( - ^) (34) 

we obtain the temperature of the black hole as 

TH = TH{l + Y.(i^^^'' (35) 

i 

where 

is the standard semiclassical Hawking temperature of the black hole and other terms are the 
corrections due to the quantum effect. Using this expression and knowing the metric coefficients 
/(r) and g{r) one can easily find out the temperature of the corresponding black hole. The same 
result was also obtained in [12] for scalar particle tunneling. 
For the Schwarzschild black hole the metric coefficients are 

/(r)=5(r) = (l-^); = 2M. (37) 
Therefore using (|35p and (j36p it is easy to write the corrected Hawking temperature: 

Th = ^{l + y H—Y^. (38) 

Now use of the Gibbs form of first law of thermodynamics gives the corrected form of the 
Bekenstein-Hawking entropy: 

[ dM 4ttAP ^ , 47rn/?2 , , , 
bh ~ J "j^ ~ — ^ °'^P^ In M h higher order terms m n 

TTrjr ^ 167r^/32 
= — h Svrpi In rjj h higher order terms in n (39) 

The area of the event horizon is 

A = Airrfj (40) 



so that 

Ah ' ' A 



S,^ = 4+4.Aln^-!^ + (41) 



It is noted that the first term is the usual semiclassical contribution to the area law ^j^jj = ^ 
[MIES]. The other terms are the quantum corrections. Now it is possible to express the quantum 
corrections in terms of S^^i by eliminating A: 

Shh = SbU + In 5bh - + (42) 
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Interestingly the leading order correction is logarithmic in A or 5'gjj which was found earlier 
in [13\ [Tl] by field theory calculations and later in [15|, [20] by quantum geometry method. The 
higher order corrections involve inverse powers of A or •S'gjj. 

To determine the value of the coefficients /?i,/32 etc we will adopt the following steps. The 
point is that nonzero values for these coefficients are related to quantum corrections (loop effects). 
Such corrections, in a field theoretical approach, are manifested by the presence of anomalies. 
Now it is a well known fact that it is not possible to simultaneously preserve general coordinate 
(diffeomorphism) invariance and conformal invariance. Retaining general coordinate invariance, 
one finds the breakdown of conformal invariance leading to the presence of nonvanishing trace 
of the stress tensor. We now show that the coefficients appearing in (j42p are related to this trace 
anomaly. 

We begin by studying the behaviour of the action ()26p upto order under an infinitesimal 
constant scale transformation, parametrised by k, of the metric coefficients. 



9fi,^ = kg^u ^ (1 + 5k)g^ 



(43) 



Under this the metric coefficients of ([TJ change as / = kf, g = k ^g. Also, in order to preserve 
the scale invariance of the Dirac equation ([2]), the field ^ should transform as ip = k^ip. On the 

3 

other hand, has the dimension of (mass) 2 and since in our case the only mass parameter is 
the black hole mass M, the infinitesimal change of it is given by, 



M = k^M ~ (1 + -5k)M. 

3 



(44) 



Now from (j33p the imaginary part of the semiclassical contribution of the single particle 
action is 



Im/i 



O(OUt) 



-2u;Im 



dr 



vWJgir) 



(45) 



where u; gets identified with the energy (i.e. mass M) of a stable black hole [8]. Therefore uj 
and h transforms like (j44p and respectively under (j43p . 

Considering only the h and ti^ order terms in (j26p and using (|44p we obtain, under the scale 
transformation, 



-^"(1+2) = ^lm-fi(out) + ^^Im-^2(out) 



Therefore 



1 , flfh\^ f 
+ ^jl"^^0(out) 

--.2 



^^'/52Vi ^ l.MT T 



^(1+2) + + 



VM2 



M4 



-5A;Im/o(out) 



(46) 



5In 



(1+2) 



leading to. 



5L 



(1+2) - -^(1+2) 



VM2 



+ 



^J-<5Mm/o(out) 



(1+2) ^ ( flPl 



5k 



Vm2 



+ 



0- 

/ 3 



^j:^im^O(out) 



(47) 



(48) 
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Now use of the definition of the energy-momentum tensor and (j48p yields, 

^ /■,4 2,5/(1+2) (hpi n^02\2 

ImJ d x^gTH = = [m2+ J^j^^'^^Hont) (49) 

Thus, in the presence of a trace anomaly, the action is not invariant under the scale transfor- 
mation. Since for the Schwarzschild black hole /(r) and g{r) are given by (j37p . from (j45p we 
obtain Im5^°"*^ = -A-kloM. Substituting this in (09]) we obtain for a; = M as, 

+ ^ = -|^Im 1 d^x^Ti: (50) 

where is calculated upto two loops. Starting from the action ()26p and following the identical 
steps as above, a similer relation among all /?'s with the right hand side of ()50p can be established. 
In this case is due to all loop expansions. 

Since the higher loop calculations to get T^j^u (from which is obtained) is very much 
complicated, usually in literature ^28j only one loop calculation for T^i, is discussed. Thus, 
comparing only the order on both sides of (j50p . we obtain, 

/?i = -|-Im J d'^xV^TH^'^ (51) 

This relation clearly shows that /3i is connected to the trace anomaly. A similar relation is given 
in |26j where it has been shown that the coefficient /?i is related to trace anomaly for the scalar 
particle tunneling. The only difference is the factor before the integration. This agrees well with 
the earlier conclusion [27\ [T3] where using conformal field theory technique it was shown that 
this Pi is related to trace anomaly and is given by, 

^1 = - - ^0 - + + - 212A^2) (52) 

3607r V 4 2 4 2 / 

'iVs' denotes the number of fields with spin 's'. In our case Ni = 1 and Nq = Ni = Ns = N2 = 0. 

2 2 

To conclude, we have successfully extended our approach [12] of scalar particle tunneling 
beyond semiclassical approximation to the model of fermion tunneling. We have considered 
all orders in the single particle action for fermion tunneling through the event horizon of the 
black hole. We showed that higher order correction terms of the action are proportional to 
the semiclassical contribution. A similar result was shown earlier in [12j for the scalar particle 
tunneling. By dimensional argument and principle of "detailed balance" the same form of the 
modified Hawking temperature, as in the scalar case, was recovered. The logarithmic and inverse 
powers of area corrections to the Bekenstein-Hawking area law were also reproduced. Finally, 
we showed that the coefficient of the logarithmic term of entropy is related to trace anomaly. 
However, the prefactor appearing in this term is different from the scalar particle example, a 
result that is supported by earlier works [271 [13] . 

Here we have only told about Discussions on other coefficients can also be given from 
(j50p . but since no information about T^i, due to multi-loops is available in the literature, we 
cannot say anything about them at this moment. 

Acknowledgment: 

I wish to thank Prof. Rabin Banerjee for suggesting this investigation and constant encourage- 
ment. 



8 



References 

[1] M.K.Parikh and F.Wilczek, Phys. Rev. Lett. 85, 50 42 (2000) |arXiv:h ep-th/ 9907001] . 
M.K.Parikh, Int. J. Mod. Phys. D 13, 2351 (2004) [arXiv:hep-th70405 160 ] . 

[2] K.Srinivasan and T.Padmanabhan, Phys. Rev. D 60, 024007 (1999) |arXiv:gr-qc/9812028] . 
S. Shankaranarayanan, K. Srinivasan and T. Padmanabhan, Mod. Phys. Lett. A 16, 571 

(2001) [arXiv:gr-qc/0 007022|. 

S. Shankaranarayanan, T.Padmanabhan and K.Srinivasan, Class. Quantum Grav. 19, 2671 

(2002) |arXiv:gr-qc/0010042|. 

S. Shankaranarayanan, Phys. Rev. D 67, 084026 (2003) |arXiv:gr-qc/0301090 1. 

[3] V.A.Berezin, A.M.Boyarsky and A.Yu.Neronov, Gravitation and Cosmology 5, 16 (1999) 
|arXiv:g r-qc/0 605099| . 

[4] M.Ar zano, A.J.M.M edved and E.C.Vagenas, JHEP 0509, 037 (2005) 
|arXi v:hep-th/050 5266) . 

A. J.M.Medved and E.C.Vagenas, Mod. Phys. Lett. A 20, 2449 (2005) 
tarXiv:gr-qc/0504113] . 

[5] Qing-Quan Jiang, Shuang-Qing Wu and Xu Cai, Phys. Rev. D 73 064003 (2006) 
|arXiv:hep-th/0512351|. 

Yapeng Hu, Jingyi Zhang and Zheng Zhao, Mod. Phys. Lett. A 21 2143 (2006) 
|arXiv:gr-qc/0611026|. 

Zhibo Xu and Bin Chen, Phys. Rev. D 75 024041 (2007) (arXiv:hep-th/0612261 ||. 
Cheng-Zhou Liu and Jian-Yang Zhu, "Hawking radiation as tunnehng from Gravity's 
rainbow" , |arXiv:gr-qc/0703055 1. 

Qing-Quan Jiang and Shuang-Qing Wu, Phys. Lett. B 635, 151 (2006) 
|arXiv:hep-th/0511123] . 

Yapeng Hu, Jingyi Zhang and Zheng Zhao, Int. J. Mod. Phys. D 16, 847 (2007) 
[arXiv:gr-qc/0611085|. 

G.E.Volovik, "On de Sitter radiation via quantum tunnehng", j arXiv:08Q3j_336 j] . 

Jingyi Zhang, "Black hole quantum tunnelling and black hole entropy correction" 

|arXiv:0806.2441j . 

[6] E.T.Akhmedov, V.Akhmedova and D. Singleton, Phys. Lett. B 642, 124 (2006) 
|arXiv:hep-th /0608098| . 

E.T.Akhmedov, V.Akhmedova, D. Singleton and T.Pihing, Int. J. Mod. Phys. A 22, 1705 
(2007) |arXiv :hep-th /0605137|. 

T.Pilling, Phys^Tett. B 660, 402 (2008) ^arXiv: 0709716241 . 

T.K.Nakamura, "Factor two discrepancy of Hawking radiation temperature", 
[arXiv:0706.2916j. 

B. D.Chowdhury, Pramana 70, 593 (2008) [arXiv:hep-th/0605197 |. 

[7] M.Angheben, M.Nadalini, L.Vanzo and S.Zerbini, JHEP 0505, 014 (2005) 
|arXiv:hep-th/0503081|. 

R.Kerner and R.B.Mann, Phys. Rev. D 73, 104010 (2006) iarXiv:gr-qc/(M30l9] . 
P.Mitra, Phys. Lett. B 648, 240 (2007) |arXiv:hep-th/0611265] . 

[8] R.Banerjee and B.R.Majhi, Phys. Lett. B 662, 62 (2008) [arXiv:0801.0200] . 

[9] R.Banerjee, B.R.Majhi and S.Samanta, Phys. Rev. D 77 124035 (2008) [arXiv: 080 1.3583] . 



9 



[10] R.Banerjee, B.R.Majhi and S.K.Modak, "Area Law in Noncommutative Schwarzschild 
Black Hole", [arXiv:0802.2i76]. 



[11] R.Kerner and R.B.Mann, Class. Quant. Grav. 25, 095014 (2008) |arXiv:0710.06T2] . 
Ran Li and Ji-Rong Re, Phys. Lett. B 661 370 (2008) [arXTv:0802.3954]. 
R.D.Criscienzo and L.Vanzo, "Fermion Tunneling from Dynamical Horizons", 
[arXiv:0803.0435]. 

Ran Li and Ji-Rong Ren, Class. Quant. Crav. 25 125016 (2008) |arXiv:0803.1410j . 
R.Kerner and R.B.Mann, Phys. Lett. B 665 277 (2008) [arXiv:0803.2246]. 
De-You Chen, Qing-Quan Jiang, Shu-Zheng Yang and Xiao-Tao Zu, "Fermions tunneling 
from the charged dilatonic black holes", [arXiv:0803.3248j. 

D.Y.Chen, Q.Q.Jiang and X.T.Zu, Phys. Lett. B 665 106 (2008) iarXiv:0804.0131j . 
Shiwei Zhou, Wenbiao Liu, Phys. Rev. D 77 104021 (2008). 
Q.Q.Jiang, Phys.Rev. D 78 044009, (2008). 

[12] R.Banerjee and B.R.Majhi, JHEP 0806 095 (2008) [arXiv:080 5.2220]. 

[13] D.V.Fursaev, Phys. Rev. D 51, R5352 (1995) [arXiv:hep-th/9412161] . 

[14] R.B.Mann and S.N.Solodukhin, Nucl. Phys. B 523, 293 (1998) |arXiv:hep-th/97 09064] . 

[15] R.K.Kaul and P.Majumdar, Phys. Rev. Lett. 84, 5255 (2000) [arXiv:gr-qc/000204 0]. 

[16] T.R.Govindarajan, R.K.Kaul and V.Suneeta, Class. Quantum Grav. 18, 2877 (2001) 
| arXiv:gr-qc/0104010 | . 

[17] S.Das, P.Majumdar and R.K.Bhaduri, Class. Quantum Grav. 19, 2355 (2002) 
| arXiv:hep-th/0111001] . 

[18] S.S.More, Class. Quantum Grav. 22, 4129 (2005) [gr-qc/041007l] . 

[19] S.Mukherji and S.S.Pal, JHEP 0205, 026 (2002) |arXiv:hep-th/0205164] . 

[20] A.Ghosh and P.Mitra, Phys. Lett. B 616, 114 (2005) |arXiv:gr-qc/0411035]. 

[21] For a review and a complete list of papers on logarithmic corrections, see D.N.Page, New 
Journal of Phys. 7, 203 (2005) |arXiv:hep-th/0409024] . 

[22] S.K.Modak, "Corrected entropy of BTZ black hole in tunneling approach", 
[ arXiv:0807.0 959j. 

[23] E.T.Akhmedov, T. Pilling and D. Singleton, "Subtleties in the quasi-classical calculation of 
Hawking radiation", jarXiv:0805 .2653] . 

[24] J.D.Bekenstein, PhD Thesis Princeton University, Princeton, NJ (1972). 
J.D.Bekenstein, Lett. Nuovo Cimento 4 737 (1972). 
J.D.Bekenstein, Phys. Rev. D 7, 2333 (1973). 
J.D.Bekenstein, Phys. Rev. D 9 3292 (1974). 

[25] J.M.Bardeen, B.Carter and S.W.Hawking, Commun. Math. Phys. 31, 161 (1973). 

[26] R.Banerjee and B.R.Majhi, "Quantum Tunneling and Trace Anomaly", [ arXiv:0808.3 688] . 

[27] S.M.Christensen and M.J. Duff, Phys. Lett. B 76 571 (1978). 

[28] B.S.DeWitt, Phys. Rep. 19, 295 (1975). 



10 



